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Abstract. The lattice potential Korteweg-de Vries equation (LKdV) is a par- 
tial difference equation in two independent variables, which possesses many 
properties that are analogous to those of the celebrated Korteweg-de Vries 
equation. These include discrete soliton solutions, Backlund transformations 
and an associated linear problem, called a Lax pair, for which it provides the 
compatibility condition. In this paper, we solve the initial value problem for 
the LKdV equation through a discrete implementation of the inverse scatter- 
ing transform method applied to the Lax pair. The initial value used for the 
LKdV equation is assumed to be real and decaying to zero as the absolute 
value of the discrete spatial variable approaches large values. An interesting 
feature of our approach is the solution of a discrete Gel'fand-Levitan equation. 
Moreover, we provide a complete characterization of rcflcctionless potentials 
and show that this leads to the Cauchy matrix form of A''-soliton solutions. 



1. Introduction 
The lattice potential KdV equation 

Qpq{x, X, X, x) — {x — x){x - x) - + = (1.1) 

where x — x{m, n), x — x{'m + 1, n), x — x(m, n + 1), x ~ x{m + 1, n + 1) and p and 
q are complex lattice parameters, is an example of a nonlinear integrable lattice 
equation. It appears as the permutability condition for Backlund transformations 
of the KdV partial differential equation (see e.g. [5]), transforms to the potential 
KdV partial differential equation under a particular continuum limit (see e.g. |14j). 
and has now been studied as an integrable lattice equation in its own right, (jl.ip 
possesses a 3D consistency property, a notion which has been studied in |15| . |12] . 
[4] and [1]. In [12] this property was shown to be equivalent to the existence 
of a Lax pair. In 2003 (jl.l[) appeared as (HI) in the exhaustive list of integrable 
lattice equations in [I] , as a representative of a particular family of integrable lattice 
equations. One-soliton and two-soliton solutions are derived in [2] by starting with 
a linear seed solution of (|1.1[) and applying the Backlund transformation of the 
potential KdV partial differential equation. Furthermore in [16] and [13] the authors 
give a determinant form for an A^-soliton solution. In this paper we solve the initial 
value problem for (jl.ip . for initial profiles satisfying 

-f CJO 



|a;™+2,o - a;,„,o - 2p|(l + |m|) < oo, (1.2) 



m— — oo 



Xm+2,0 - Xmfi > (1.3) 

along n = 0. The m and n variables are related to the " and ' shifts hy x = Xrn+i,n 
X = Xm,n+i- While it may appear that equation (jl.ip only depends on the single 
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parameter — q^ , it is natural to maintain the dependence on both p and q because 
each is associated with a different Backlund transformation, and one is therefore 
able to isolate the effect of each individual transformation. The initial conditions 
(|1.2p and (|1.3p arc posed along the axis n — Q, which is analogous to the continuous 
case, in which the initial condition is posed on i = 0. This is the reason why the 
parameter q associated with the n-dircction does not appear here. The summability 
condition (jl.2|) is a direct analogue of the intcgrability condition placed on the initial 
condition in the continuous case, and imposes the condition that the asymptotic 
behaviour of Xm+2fi — Xm,o — 2p decay faster than as |m| — > oo. Analogously 
to the continuous case, one finds for soliton solutions that this decay is exponential. 

1.1. Background. To our knowledge the first studies into the discrete inverse 
scattering transform date back to Case and Kac [6] and Case [7]. These authors 
considered a direct discretisation of the time-independent Schrodinger equation and 
were led to an eigenvalue problem A(p = where A is a tridiagonal matrix. This 
was solved as an initial value problem for the half- line n > 0. The inverse problem 
was posed on the unit circle in the complex plane of the spectral variable, and 
using the orthogonal polynomials that arose from the spectral distribution of A the 
solution was obtained by deriving a discrete Gel'fand-Levitan integral equation. 
Flaschka showed how this procedure could be applied to solutions of the Toda 
lattice. The author considered linear difference equations in which the coefficients 
depended on the Hamiltonian of the lattice. These coefficients were assumed to 
depend smoothly on time and thus the discrete spectral data evolved according to 
a continuous evolution equation. The eigenvalues were shown to be constants of 
motion for all time, and the solution of the inverse problem was given as a solution 
of a discrete Gel'fand Levitan integral equation. 

More recently Boiti et al [5] considered an "exact"Q discretisation of the Schrodinger 
equation. They were led to a different discrete problem to , [7] and [9] , one which 
had been studied earlier by Shabat jl8], and coupled this with a two-parameter 
differential-difference time evolution equation, reminiscent of that given in the con- 
tinuous inverse scattering transform in [2]. Rather than following the Gel'fand- 
Levitan method for the inverse problem, the solution was given in terms of the 
time-dependent Jost solutions in the expansion of the spectral variable. These au- 
thors also considered a discretization of time, and postulated a similar evolution 
equation as in the continuous-time case. They found that this particular discretiza- 
tion led to a higher-order version of (|l.ip . 

This spectral problem in [5] was later considered by Shabat [l^ as a dual problem 
to the spectral problem for the continuous Schrodinger equation. Using this duality 
property the author was able to use results of the continuous inverse scattering 
transform to obtain qualitative estimates for the discrete spectral problem. Levi 
and Petrera |10| obtained this spectral problem as one of the Lax equations for 
(|l.ip and used it to solve the inverse scattering problem for (|l.ip . Their discrete 
"time" evolution was obtained using the second Lax equation and was different to 
that presented in [5], however the inverse problem was again given in terms of the 
Jost solutions and thus was a generalisation of that given in j5j . 

In 2002 Ruijsenaars [T7] considered the same linear discrete problem as [5], cou- 
pled with a different parametrisation of the associated differential-difference time 



Exact" meaning that the discretization arises directly from applying Darboux transforma- 
tions to the Schrodinger equation 
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evolution equation. This equation depended on two independent variables x and t 
and contained iterates in x and derivatives in t. Ruijsenaars assumed analyticity 
in X, and considered the discrete problem as an eigenvalue problem for an analytic 
difference operator. This was preceded by a study of a class of rcflectionless ana- 
lytic difference operators that yielded soliton solutions. The inverse problem relies 
on a Hilbert transform, and the soliton solutions arc shown to converge to regular 
KdV solitons under suitable scaling limits. 

In contrast the work presented here deals with the solution of the initial-value 
problem for (|l.ip as a function of two discrete independent variables. Thus both 
the linear problem and the scattering data evolve according to difference equations. 
We do assume continuity in the lattice parameters which allows the inverse problem 
to be solved by a Riemann-Hilbert approach, and we show that solutions of 
are given in terms of solutions to a discrete Gel'fand-Levitan integral equation. 



1.2. Outline of Results. In this paper we solve the initial value problem for the 
LKdV (|l.ip rigorously through the inverse scattering transform method. Following 
|12j we obtain a Lax pair for and derive the governing linear equation which 
agrees with that in [10]. In Sections (jSj, (HJ and ([5)) the direct scattering procedure 
is carried out, in which the analyticity properties of the Jost functions in the plane 
of the spectral variable are proved rigorously, obtaining more precise bounds than 
those given in [5] . We also give a sufficient condition for the poles of the transmission 
coefficient to be simple. In [5] and |10| this property was assumed without proof, 
however this is not true in general for an arbitrary potential as is seen in the given 
counterexample. In Section ([6|) the discrete "time" evolution of the transmission 
and reflection coefficients is then derived, and the inverse problem is treated in 
Section ([7]). Rather than give the solution in terms of the "time" -dependent Jost 
solutions, as was done in [5] and [10], we emulate the Gel'fand-levitan procedure 
for the continuous case and derive a discrete Gel'fand-Levitan integral equation 

ni 

K{m,L) + B{L) + K{m,r){B{r-m + L) + B{r-m + L-l)) = 0, 

r— — oo 

where B is dependent on the scattering data. The solution K(m^L) of this hnear 
equation is then related to the solution of by 



1 + K{m + 2,m + 2) 
1 + K{m + l,m + 1) 



In Section (jSJ) we consider one- and two-soliton examples and finally in Section ^ 
we show that any rcflectionless potential that satisfies the required summability 
and positivity conditions gives rise to an iV-soliton solution identical to that found 
by applying Backlund transformations [13] . 



2. Lax pair 

Nijhoff 's method for obtaining a Lax pair for the Adler system [12] is amenable 
to (jl.ip . and relies on the multidimensional consistency of the equation. If a; = 
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x(rn, n,l + 1) denotes a shift in the lattice with parameter r then we have 

Qpq{x, X, X, x) = (2.1a) 







(2.1b) 
(2.1c) 



along with the group of square symmetries Qpq{x^x^x^x) — Qpq{x^x^x^x) — 
Q qp(^x ^ x^ x^ x^ . Since pTTj) is fractional linear in each variable we may solve (|2.1b|) 
and (|2.1cp for 5; and x respectively 



— x(x — x) + 
X = z 



X — X 
x{x — x) + q'^—r'^ 



(2.2) 
(2.3) 



The fact that (j2.1bp and p.lcp arc both discrete Riccati equations for x suggests 
the separation x = f/g, so x ^ f/g and x = f/g. By defining 4> ^ 
()2.2p and p.3p can be written in matrix form as 

(f) = KlLcf) 



4> = K2AI(f) 



where 



L = 



M 



-X XX 
-1 



■P 

X 



-X XX + — 
-1 X 



then 

(2.4) 
(2.5) 

(2.6) 
(2.7) 



and Ki,K2 are as yet undetermined separation functions. The determinantal con- 
dition 112] on these functions is then 



(KiK2)"'(detL)(dctM) = (KiK2)"'(det M)(det L) 



(2., 



Since detL 



P 



and detM 



this allows for scalar (w.l.o.g. unity) 



values of the separation functions. One then finds that 

^ - = Qpq{x,X,X,x) 



X + xj 

-1 







(2.9) 



so that = necessarily implies x solves (jl.ip . 

The two systems of first-order difference equations p.4p and p.Sp for / and g 
give rise to the second order difference equations for g 

g — {x — x)g + {p^ — r'^)g = 

g-{x~ x)g + {q^ - r^)g = 0, 



(2.10) 
(2.11) 



along with f ~ xg ~ g and f ~ xg — g. Equations ()2.10p and p. lip serve as a 
Lax pair for the inverse scattering transform, since from (|2.9p we have g — g = 
^Qpq{x,x,x,x)g. The lattice parameter r acts as the spectral variable and the 
difference x — x acts as a potential. The difference equations in / are auxiliary to 
the problem and are not required to be solved explicitly. 
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3. The Direct Scattering Problem 



In this section wc carry out the direct scattering procedure for equation (|2.10p . 
where we assume the solution x of p.ip to be reaL We first define r := iz, so that 
z now acts as the spectral variable, and alter the notation such that g = gim, n; z), 
g = g{m+ l,n;z), g = g{m,n+ l;z) and x = Xm,n, i = Xm+i,n, x = Xm,n+i- 

3.1. Motivation. One-soliton solutions to (|l.ip are given in [3] and [13]. In [3] the 

solution was obtained by taking a Backliind transformation of a linear seed solution 
and then using (jl.ip to obtain the m and n dependence, while in |13j the authors 
began with a Cauchy matrix structure and showed that it solved the homogeneous 
version of (jl.ip . The one-soliton solution is given by 



pm + qji + C + 



2k 



where 



Pp = 



p + k 



Pq 



q + k 



(3.1) 



(3.2) 



p — k q ~ k 

and C is constant. Thus for any n the leading order behaviour of this solution as 
m — > ±cxD is 

Xm,n ^ pm + qn + const, as, m ^ ±oo, 

where all other lower-order terms vanish exponentially. One finds the same result 
for the two-soliton solution given in Section ([5]) . We therefore assume 



where u is a real- valued function satisfying u — > as m 



(3.3) 

±oo, independently of 



3.2. Initial conditions. As n is arbitrary in the direct scattering problem we set 
n = and define gijn; z) = g{m, 0; z) and Um = Um.o- Since (jl.ip is invariant under 
the maps p H- —p and (/ i— > —q we set p > and g > without loss of generality. 
We also assume that 

2p + Um > for aU m, (3.4) 

which is satisfied by p.ip provided we have A > and p > k > 0, and by the 
two-soliton solution (with a similar restriction on parameters) given in Section ([S]) . 
This assumption is sufRcienlQ to prove that all discrete eigenvalues are simple (which 
was assumed without proof in |S] and jlOp. and as a consequence implies that all 
discrete eigenvalues are purely imaginary and lie within the interval (0, ip). 

With the above assumptions the direct scattering problem is entirely governed 
by the second-order difference equation 

g{m + 2;z)- {2p + u,„+i).g(m + 1; z) + {p^ + z^)g{m; z) = 0. (3.5) 



^but not necessary. If u,n = c{5m,o + 5m, i) then the case c = — p(l + \/2) gives 

z(z + ip) 

where the zeroes of a(z), which is defined by 1 13.71 1. are the discrete eigenvalues. For c < — p(l-l-v^) 
however, the discrete eigenvalues arc simple. 
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3.3. Jost Solutions. 

Definition 3.1. The Jost solutions {(p,'ip) and (ipjil^) to p.Sp arc defined by tlie 
boundary conditions 

ip ^ {p- iz)"" /o r \ 

_ ; , . N„, as TO — >■ —oo (3.Da) 

If (p + iz)"^ ^ ' 

%■ /^~'~* Nm as TO — +CXD, (3.6b) 

■0 ~ (p - IZ)™ ^ ' 

where (j) ± iz)™ solve p.Sp as u — > 0. 



3.4. Spectral Data. Equation (j3.5p is invariant under the transformation z — z, 
and by considering how the boundary conditions of the Jost solutions change under 
this mapping, uniqueness of the boundary value problem (see e.g. implies that 
^(to; z) = ip(m,\ ~z) and ■(/'(to; z) = ?/'(to; — z). Since the general solution involves 
two linearly independent solutions we may write 

■0 = cup + h(p^ (3.7) 

where a = a{z) and b = b(z), and thus ip = a(—z)ip+b{—z)if. If z is real then taking 
the complex conjugate of p.Sp reveals Tp = ip* and ^ = if)* , and so -0 = a(/3* + 6(y9 
on = 0. If we define the general solution to (|3.5p by the boundary conditions 

- - 1 + R{z) ( ™ as TO ^ -oo (3.8a) 



[p + iz)'' 
.g(TO; z) 



W^{gi,g2]z) = ^ ^ „^ [gi(TO;z)g2(TO + l;z) - 5i(to + 1; z).g2(TO; z)] . (3.9) 



r(z) as TO +0O, (3.8b) 

where R and T are the transmission and reflection coefficients respectively, it follows 
that g = lp + Hp and g = Ti' = TaJp + Tbip. Thus i? = ^ and T = i. 

Definition 3.2. Let gi and 172 be two solutions to p.Sp . The discrete Wronskian 
Wm{gi, g2', z) of gi(m; z) and g2{m\z) is then defined to be 

1 

[p-^ + z'^ 

Lemma 3.1. Suppose gi and g2 are linearly independent nonzero solutions of p.Sp . 
Then Wmigi, g2] z) is both nonzero and independent ofm. 

Proof. Since gi and 52 both satisfy (j3.5p we have 

51 (to + 2;z).g2(TO + l;z) + (p^ + z'^)g2{m + 1; z)gi{m;z) 
={2p + Ura+i)g2{m + 1; z)gi{m + 1; z) 

=g2(TO + 2; z).gi(?7i + 1; z) + [p^ + z'^)gi{m + 1; z)g2{m;z). 

Thus 

gi{m;z)g2{m + 1; z) - gi(TO + 1; z)g2{m; z) = Co(p^ + z^)™ 

where Co is constant as required. If Wm{gi,g2]z) = then clearly gi and 172 are 
linearly dependent. □ 

Proposition 3.1. For 5z = the functions a and b defined by (|3.7p satisfy 

\a\^-\b\^ = l. (3.10) 



AN INVERSE SCATTERING TRANSFORM FOR THE LATTICE POTENTIAL KDV EQUATION 



Proof. The linearity and anti-symmetry of the Wronskian imply that since ip = 
atp* + bip along 5z = 0, 

WUr,^;z) ^ (|a(z)|2 - \biz)\^)W^{cp,ip*;z). 

As Wmi^'* 1 'ip] z) and Wm{f,'P*',z) are independent of m, they are equal to their 
boundary values at to — > — oo and m — +00 respectively, which is 2iz in both 
cases. □ 



4. Analytic properties of the Jost solutions 
For the subsequent analysis it is convenient to make the following definition 
Definition 4.1. 

(p{m;z) 



X(m;z) := 
T(m; z) := 



[p — izY 
{p + izY 



(4.1) 
(4.2) 



Then x ^ 1 and T — > 1 as to — >■ —00 and m +00 respectively, x ^-nd T are 
defined by x('^; z) ~ x(™; ^-nd T(m; z) ~ T(m; ~z). 

Lemma 4.1. For > 0, z ^ the functions x o.'^^d T satisfy the following 
summation equations: 



1 



m — 1 



x(to;z) = 1 + — 



+ CXD 



p — iz 

p — iz 
p + iz 



m-j 



- 1 



1 



UjT{j-z). 



(4.3) 
(4.4) 



J— m+l 

Proof. Equation (|3.5p for x(™; ^) gives 

(p - iz)x{in + 2;z)~ 2px{m + 1; z) + {p + iz)x{m; z) = Um+ixim + 1; z) 
which can be summed from an arbitrary Mq to to — 1 > Mq to give 

m — 1 

{p-iz)[x{m+l; z)-xiMo+l; z)]-{p+iz)[x{m; z)-x{Mo; z)] = ^ Uj+ix{j+l; z). 

j=Mo 

Letting Mq — > —00 and incorporating the boundary behaviour of x("^; z) gives 



p — iz 
p + iz 



x{m + l;z)- x(m; z) 



-2iz 



p + iz p + iz 



^ m 



] = ~co 



/ _• \™ 

We now multiply this equation by the summing factor (2^1 and sum from an 
arbitrary Mi to to — 1 > Mi. By letting Mi ~oo, incorporating the boundary 
behaviour of x and noting that < 1 for 3z > 0, this yields 



X(m;z) = 1 + — V 



p + iz 



p + 12 ' \p — IZ 

I— — 00 



m-l I 
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Changing the order of summation then gives 

m — 1 



X(m;z) = l + — i— J2 ^3XU;z)Y, 



m~l / s m — l 

f P + IZ 



1 



— F 

2iz /-^ 



J = -QC 

m—l 



p + iz 
p — iz 



\p — iz ^ 

1 "ixO»- 



j = -oo 

The proof of (j4.4p follows by a similar argument. 



□ 



Lemma 4.2. For z = the Jost solutions x(m; 0) and T(m; 0) satisfy the following 
summation equations: 



m — 1 



X(m; 0) = 1 + - V (m - j)ujxU'^ 0) 
T(m;0) = l + - ^ (.7-mKT(j;0). 



Proof. 



(4.5) 

(4.6) 

□ 



Proposition 4.1. For SJz > 0, z 7^ the summation equations (|4.3p onrf (|4.4I 
/ia?;e i/ie following series solutions: 



+00 



fe=0 

^ Jk{m;z) 
T{m;z)^}_^ -k 



fe=0 



where 

i?o(m;z) = 1, Hk+i{m;z) 



7n — 1 



2i ^ 

Jo{m; z) = 1, Jfc+i(m; X! 



] = -oo 
+00 



j—m+l 



p + iz 

p — iz 

p ~ iz 
p + iz 



- 1 



(4.7) 
(4.8) 

u,Hk{j;z) (4.9) 
u,Jk{j\z). (4.10) 



Proof. Inserting (|4.7p into the summation equation (|4.3p gives 



m — l 



X{m-z)^l + — J2 



] = ~oo 



p + iz 
p — iz 

V ^ ( 1 V 



k=0 



- 1 



p + iz 
p — iz 



u 



+00 

E 

fe=0 

m-j 



Hk{j;z) 



UjHk{j;z) 



1+E 



Hk+ijm; z) 



fe=0 



as required. The proof of (|4.8p is similar. 



□ 
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Proposition 4.2. For z ^ the summation equations ()4.5p and ()4.6p have the 
following series solutions: 

xK0)==E%^ (4-11) 

fe=o ^ 

T^KO) = E%^ (4.12) 

fe=0 ^ 

where 

m — 1 

i/0(m) - 1, i/^i(m) - E ("^ - jhjH'kU) (4.13) 

J°(m) = 1, Jfc+i(m; z) = ^ (^ " ^>jJkij)- (4-14) 

Proo/. □ 

The following theorems describe the analyticity properties of the Jost solutions 
in the z-plane. These results closely mirror those obtained for the direct scattering 
of the continuous Schrodinger equation, which can be found in the detailed analysis 
given in [5]. 

Theorem 4.1. Assume 

E 1^^,1(1 + |j I) <oo, (4.15) 

j=-oo 

Then for > 

|x(m;z) - 1| < Ci for z 7^ (4.16) 
|x(m; z) - 1| < C2 (1 + max{m, 0}) , (4.17) 

(4.18) 

where Ci and C2 are constant. For all m, xl'^i; z) is analytic in > and con- 
tinuous in 3z > 0. For all 3z > (|4.7p converges absolutely in m (and uniformly 
zfzj^O). 

Proof. In the following all symbols Ni refer to constants whose precise values are 
not required, but are used to obtain the required results. We first prove (|4.16p . For 
Sz > we have < 1 and so the recursion relation (|4.9[) for x(™; z) can be 

upper-bounded by 

rn — 1 

\Hk+i{m;z)\< J2 \uj\\HkU;z)\. (4.19) 
\Hkim;z)\ < ^ ^, ' 

where 



j = -oo 

Lemma 4.3. 



j = -oo 
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Proof. Clearly this holds for fc = 0. Equation (|4.19p then implies 



\Hk+i{m;z)\< 



kl 



] = -oo 

m — 1 



^ E iP{j)-P(j-i))P{j~iY 



Summing by parts then gives 



3 = -ca 



— OQ j — ~QO 

rn-l /k-1 \ 

p{m - 1)^+1 - E pumj) - pu - 1)) X E - !)'■ 



] = -co 

+ CXD 



□ 



<p(m-if+i-fc E (n?) - n? - i))n? - 1)' 

since P{j — I) < P{j) for all j. This completes the inductive step. 
By (|4.7p and Lemma (|4.3p we have 

which gives (|4.16p . Thus for any z satisfying 5z > 0, z 7^ the series solution for 
x(m; z) converges absolutely and uniformly in m. An alternative upper bound for 
(1431) is 



|x(m;z)| < 1 + 



1 



\p-iz\ 



E 



E 



p — IZ 
p + iz 



where a = max{l, - }. Thus a majorant for either (|4.3p or (|4.5p is 



7n — 1 



|x(m;z)| < 1 + cr E - J)l%-Ilx(j;^) 



j = -oo 



Therefore for Sjz > 



where 



+00 



|x(m;z)| <EHfe(m) 



(4.20) 



(4.21) 



fc=0 



Hfc+i(TO)=cr E ("^ - J>jHfe(j)- 

j=-oo 

Following Lemma (|4.3p one can then show that 

(7^Q{m — 1, TO - 1)*^ 



|Hfc(m)| < 



fc! 
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where Q{m,L) =■ X]j'=-oo("^ ^ J + 1)1% I- Thus by ()4.2ip we have 



\x{m; z) -\\<2_^ — — - ^Q{m^ m)e^^ ' > . 



k=l 



k\ 



For m < equations and give 



|x(m;z) - 1| < ae'^«("^°) 



m— 1 

j = -oo 

-1 



j = -oo 



——00 

E i-^-) 



3^-00 



For m > equations f|4.2Up and f|4.22p give 



m — 1 



rn — 1 



|x(m;z)| < 1 + 0- E (-J)|wjl|x(j;2;)l + WCT E Wj\\x{j\z)\ 
J— —00 J— —00 

-1 rn-l 

<i + o' E ("•?)l"jllx(j;2:)l + wo- E l"jllx(j;^)l 



j = -oo 

m—1 



<N3 + ma E l"illx(j;2)l- 

j = -oc 

By writing x(to; z) = A'3(l + m)Si(m; z) we see 

rn— 1 

|Si(m;z)| < 1 + a E (1 + IjDIu.piO'; z)|, 

j = -oo 

which can be iterated to give 



|Si(m; 2)1 < cxp cr E + 



j = -oo 



Thus for m > 0, g^, (|4:22)) and (liTM)) give 



|x(m;2;) - 1| < (7 ^ (-j)|Mj||x(j; z)| + Jntr ^ |uj||x(j;z)| 

J — — 00 J^ — 00 

771—1 

<N5+mcrN4 E + 



j = -oo 



< Ne{l+m). 



(4.22) 



(4.23) 



|x(m;z)| < A'"4(l + to). (4.24) 



(4.25) 



Combining this with the upper bound (j4.23p for m < gives (|4.17p . This estimate 
is vahd for all m and all z satisfying Sz > 0. For every m the series solution 
(|4.3p for x(™; z) converges absolutely and uniformly in z satisfying SJz > 0. Thus 
x(to; z) is continuous in z within 3z > 0. Since the iterates Hk{m\ z) are analytic 
functions of z in 3z > 0, x('^; ^) is also analytic in this region. For z = the series 
solution for x("T'; 0) converges absolutely in m and uniformly for m < mo- 
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□ 

Remark 4.1. By considering the series solution for T(m; z) one can similarly prove 
its existence and continuity in > and analyticity in 3z > 0, provided ()4.15p 
holds. One obtains similar estimates to (|4.16p and (j4.17p . 

Theorem 4.2. Assume that 

^ (l+/)|u,|<^. (4.26) 

j = -oo 

Then for '^z > 

\x'{m;z)\ < (73(1 + mmax{m, 1}), (4.27) 

where C3 is constant. For all m, x'("^; ^) exists and is continuous in z for all z in 
3z > 0. 

Proof. Again let Ni denote constants as necessary. We rewrite (|4.3p as 

X(m; z) = IH ( ) UjxU] z), 

which agrees with (|4.1ip at z and is therefore valid everywhere in 9z > 0. 
Taking an upper bound of the derivative of this equation then shows 

m — 1 m — 1 

|x'(m;z)|<a2 ^ (m - |u,||x(j; z)| + a ^ {m ~ ])\u,\\x' {r, z)\. (4.28) 

j=~oo j=-oo 

If m < then (|4Tf|) implies 

m— 1 m— 1 

j — — QO j — — oo 

If m > then f|4.17p gives 

m— 1 m— 1 rn— 1 

(m- j f\uj\\xij;z)\ <2 Y m^Wj\\xij;z)\+2 ^ f\uj\\xij; z)\ 

j=-x i=-oo j = -oo 

— 1 m— 1 ni—l 

j = -oo j = l j = -oo 

m — 1 

j=-oo 

< 7V4(l + m2). 
Therefore (|4.28p becomes 

m — 1 

|x'(to;z)| < cr^iV4(l + mmax{TO, 0}) + (T ^ (m - j)|uj||x'(j; z)| (4.29) 

j = -oo 

which can be iterated to give 

|x'(m;z)| < a^N4{l + mmax{m,0})e'"^^"'-^'"'-^l 
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Inserting this into (|4.29p then gives 

|x'('Ti; z)\ < a^Ni (1 + mniax{m, 0}) 

-1 m-l 

j — — 00 j— — QO 

ra— 1 

< A^5 (1 + m niax{m, 0}) + cftii |x'(j ; ^)|- 

For m < 

m — 1 

|x'(to;z)| < iVs + CTTO ^ \uj\\x'{j]z)\, 

j=-oo 

and so by defining x' i'"^] z) = -^5(1 + |'7i|)S2(m; z) we have 

m — 1 

|S2(to;z)| < 1 + CT ^ (1 + |j|)|uj||S2(j;z)| ^ |S2(m; z)| < A^g, 

j=-oo 

and so |x'(to;z)| < Nj{l + |m|). For m > let x'{'m\z) = iV8(l + 771^)23(771; z). 
Then 

m— 1 

|S3(m;z)| < l + a ^ {I + f)\uj\\^^{j- z)\ ^ 123(777; z)| < iVg, 

j = -oo 

and so |x'(77i;2;)| < A'g(l + 777^). This proves (j4.27p . For every 777 the summation 
equation for x'(777;z) converges absolutely and uniformly in z satisfying 5z > 0, 
and absolutely in 777. Therefore x'l"^! z) exists and is continuous for all z in 3z > 0. 

□ 

Remark 4.2. Using the summation equation for T one can similarly prove the 
existence and continuity of T'{m-. z) in 3z > provided (|4.15p holds. One obtains 
a similar estimate to (|4.27p . 

Corollary 4.1. Assuming (j4.15p holds, the Jost solutions x and T exist and are 
continuous in < 0, aTid are analytic in Qz < 0. // ()4.26p also holds, then the 
derivatives of these functions also exist and are continuous in 3z < 

Remark 4.3. The existence and continuity of the derivatives of the Jost solutions 
along = is not strictly required. Theorem (|4.2I) merely serves to illustrate 
under what restrictions on u this condition will hold. Henceforth we only assume 
that u satisfies (j4.15p . 

5. Analyticity and Asymptotic Properties of a(z) and b{z) 

We now consider properties of the functions a and b, which represent the trans- 
mission and reflection coefficients by the relations R — ^ and T — ^. In particular 
we look at their analyticity properties and asymptotic behaviour as |z| +00. We 
assume that (|4.15|) holds. 

Proposition 5.1. The functions a{z) and b{z) defined by equation (j3.7|l have the 
following properties: 

- a(z) is analytic in the region 3z > 0, and is continuous on 3z = 0, except 
possibly at z ~ 
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- b{z) is continuous on = 0, except possibly at z = 0. 
Proof. We have 

W™(<^,-0; z) = Wjn{'f,a(p + b(p; z) = a(z)W™(v3, ^; z) = 2iza{z) 
Wmii','^, ; z) = Wm{aip + bip,7p; z) = b{z)Wra{^p,W\ z) = 2izb{z). 

Since and ip are both analytic in 3z > and continuous on = 0, a{z) also 
has this property, except possibly at z = 0. The expression for b{z) however is 
only valid on 32 = 0, where </? and ■0 are both defined. Thus b{z) is continuous on 
3z = 0, except possibly at z = 0. □ 

Proposition 5.2. For z ^ Q the Junctions a{z) and b[z) can be expressed as 

a(z) = l-— u,T{j;z) (5.1) 



2iz 

j = -co 



1 / -I- ■ \ 



j = -oo 



Proof. The summation equation (j4.3p for x can be expressed as 



T(m: z) 



\ j=m+l J 

( 1 



Comparing this with T{m] z) ~ a{z) + b{z) ( J as m — oo gives the desired 
result. □ 



Proposition 5.3. 



X(m; z) - 1 + O ( i ) as |z| +oo with 5z > (5.4) 



T(m; z) - 1 + O as |z| ^ +oo with 3z > 0. (5.5) 

froo/. This is clear from (|iT7)) and □ 
Corollary 5.1. 

X(m; z) - 1 + O as |z| +oo with 5z < (5.6) 

T{m;z) r^l + O (-j as |z| +oo with Sz < 0. (5.7) 
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Corollary 5.2. 

a(z) - 1 u.j+0{ — ] as Izl -> +00 with SJz > 0, (5. 



2iz ^ ■> 

j=-oo 



j — — CO ^ ^ 



Proof. This follows directly from inserting the asymptotic behaviour (|5.5p of T into 
the expressions (|5.ip and (|5.2p for a and 6. 

□ 

Theorem 5.1. a(z) has a finite number of zeroes Zk, k = 0,1,..., N in the open 
half-plane '^z > 0, and at each Zk we have ip(m; Zk) = bk(p{m; Zk), where bk = b{zk). 
Moreover every zero is simple, lies on the imaginary axis and satisfies \zk\ < p. 

Proof. Proposition (|3.ip implies > 0, and since a(z) ^ I + O (i) as \z\ — > 00 
it follows that each Zk must be finite in magnitude. It is also true that since 
a{z) = ^Wmi'P, tp', z), Lemma p.ip implies that iplm; Zk) and ^[m; Zk) are linearly 
dependent, so we write ip{m; Zk) — bky^im; Zk), for some constant bk. This implies 

T(m; Zk) ^ bk — ) as m ^ -00, 

\p + lZkJ 

so by equation (|5.3p we have 

+00 

2izk = ^ UjT{m; Zk) 

j=~oo 

2iZk \p-izk, 
J— —00 ^ ' 

and thus bk = b{zk). 

To show the zeroes of a(z) are purely imaginary, consider 



(jf2 _j_ 2 



At z ~ Zk we have 



P - IZk 

P + izk 

P + izk 
P - izk 



— > as m — ^ ~oo 



as TO — > +00, 



and so X^jLT^oo '^(j! -^fc) exists. Equation p.Sp for v becomes 

(/ + 2^)5 (iy(m + 2; z) + i^{m; z)) = (2p + u,„+i)i/(to + 1; z), (5.10) 

and by fixing z ^ Zk, multiplying by v*{m + 1; z^) and summing each term from 
TO — —00 to +00 we obtain 

+00 

{p' + zl)i{S + S*)^ E {2p + u,+,)Hj,Zk)\', 



16 SAMUEL BUTLER, NALINI JOSHI 

where S = I^j^oo ^0 + 1' Zk)y*{j; z^). Since we assume (2p + Um+i) > for all 
m it follows that the right hand side of this equation cannot vanish, and so zf. must 
be real. This in turn implies that il){m]Zk) is a real-valued function. Hence by 
squaring (|5.10p and fixing z = we see 

(p2 + zl) [v{m + 2; Zfe) + v{m; Zk)f = [{2p + it„+i)i/(m + 1; Zk)f > 0, (5.11) 

which implies that — |zfcp > 0. 

To show that the zeroes of a{z) are simple we note that since a{z) — ^Wmi'P, ip', z) 
it follows that 

2iZka'{zk) = ^W{^,^-z)\ =W{^',^;zk)-W{,p',^;zk). (5.12) 

dZ IZ-2fc 

Taking a z-derivative of p.5p shows that ip' satisfies 

(p'{m. + 2;z) + {p^ + z'^)ip'{m; z) + 2z(p{m; z) — {2p + Um+ij'-p' {m + 1; z) 

^(p'{m + 2; z)^jj{m + 1; z) + [p^ + z'^)ip{m + 1; z)(p'{m; z) + 2z%jj{m + 1; z)Lp{m; z) 

= {2p + u„i+i)ip'{m + 1; z)ip{in + 1; z). 

Comparing this with equation p. 51) for ip shows 

(p'{m + 2; z)^{m. + 1; z) + {p^ + z'^)^jj{m + 1; z)(p'{m; z) + 2zil>{m + 1; z)(p{m,; z) 

= ipim + 2; z)Lp'{m + 1; z) + {p^ + z'^)ip'{m + 1; z)il;{m; z), 

which after dividing through by (p^ + z^)'"+^ and summing from mo < m to m 
reveals 

The same procedure can then be carried out with (p and ip interchanged, except 
choosing to sum from m to mi > m: 

m^,^;^)|„, -2z)_^ („2 + ^2y+i ■ 

j"=m 

Now fix z = Zfc. Since ip{m] Zk) ~ bkf{m; Zk) it follows that 

IP — IZl^ \ 

W{ip' , 'ipj Zk) ^ ibk I ^ as m ^ — oo 

\p+lZkJ 
— Z I \ '^Zk \ 

W{ip', ip; Zk) ^ — — ^ as m ^ +cx), 

bk \p - izk J 

and so by taking mo ~ —oo, mi = +oo, equation (|5.12p shows 

'( ^ - V'(j + l;zk)ipij;zk) 

"^''''bk.^ (p2 + z2),+i • 

From p.Sp however we have 

V'(m-|-2; 2;fc)'0(m + l; Zfc) + (p^ + Zfc)^/'(m + l; Zk)ilj{m; Zk) = (2p+w™+i)V'^(m-t-l; Zfc), 
which implies 

^ (2p + Uj-+i)^2(j + l;zfc) 

^ ~ fe^ 2(p2 + z2).+i ■ ^^-^^^ 

J — — OO ^ 
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As '(/'(™; Zk) is a real-valued function this final sum, and hence a'{zk), is nonzero as 
required. 

Finally since a{z) is analytic in 5z > 0, it must have isolated zeroes in a finite 
region along the positive imaginary axis. The only way that a{z) could have an 
infinite number of zeroes in this region is if they formed a limiting sequence which 
accumulated at z = 0. We show that this is not possible. Assume that there exists a 
sequence of zeroes of a(z), {z^.}, which all lie on the imaginary axis and accumulate 
at z = 0; limfc^.oo Zk = 0. Then at each z^ we have 

'(/'(m;zfe) 



b{zk) 



and so 



lim \b{zk)~bm 



lim 



ip{m;zk) ' 
-ip{m;zk) 



ip(m] 0) 



k~^oo k- 

since the Jost solutions are continuous at z = 
ip{m; 0) = ^(m; 0) and so tpim; 0) = (a(0) + 6(0))i^(m; 0) 

lim 6(zfe) = lim ^''^ 

k^oo k~yoo ip[m', Z)^) 







(/5(m;0) 

At z = however we have 
But 



6(0): 

a(0) 



a(0) 
0, 



which contradicts the fact that |a(0)p = 1- 
of zeroes in Sz > 0. 



6(0) 



' |6(0)p > 1. Thus a has a finite number 

□ 



6. Evolution of the Scattering Data in the ti-direction 

Equation (f2ll|) governs the behaviour of solutions in the n-direction: 

g{m, n + 2; z) - {xm,n+2 - Xm,n)gim, n + 1; z) + (g^ + z^)g{m, n; z) = 0. (6.1) 

The one- and two-soliton solutions behave like Xm,n ^ pm + qn+const. as m ^ ±oo 
independently on n, where all lower-order terms vanish exponentially. We therefore 
assume that Xm,n+2 — Xm,n — ^ 2q as m ±cxd for all n. 

Definition 6.1. The n-dependent Jost solutions (^'■"^ and '(/'*•"•' to equation (|6.1|) 
are defined to be 

(^(") = (g - iz)"(p (6.2) 
= (9 + iz)"V- (6.3) 

^'•"^ and ^ are defined by Tp'^'^^m^n; z) = ip^"^ {m,n; —z) and \m,n;z) = 
■0'^")(to, n; —z). 

The n evolution of reveals that ip satisfies 

{q'~iz)ip{m,n + 2] z) - (x,„,„+2 -2:m,r^)v("^, + 1; z) + {q + iz)(p{m, n; z) = 0. (6.4) 

Note that (|6.ip does not allow for the ^-independent boundary conditions of ip 
and ip, but these are consistent with (j6.4p . Since 'ip{m, n; z) — a{n; z)Tp[m, n; z) -|- 
6(n; z)ip{m, n; z) holds for all n, we have 

V'^"^"^,^;^) = a(n;z)^(")(m,n;z) + 6(n;z) f^-t^^ ^Jj'-''\m,n; z), (6.5) 

\q - izj 
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which agrees with p.7|) on n = if we adopt the notation a(z) = a(0; z) and 
b{z) = 6(0; z). 

Theorem 6.1. The function a{n; z) is independent of n: 

a{n; z) = a{0; z) = a{z). (6.6) 
The function b{n; z) is given by 

Proof. In the hmit m — — oo we have (p^''^'^ {m,n; z) ^ {p ~ iz)™{q — iz)^ and 
^'•"•'(m, n; z) ~ (p + izY'^{q + iz)", and so in this hmit i^'"' and ^^^"^ are hnearly 
independent solutions of (|6.ip . Thus we may write 

?/'^"^(m, n; z) ^ Ci(m; z)^*-"''(m, n; z) + C2(?ti; z)93''"^(m, n; z) as m — > — oo, 

where Ci and C2 are independent of n. Comparing this with (j6.5p gives the desired 
resuh. □ 

7. Inverse Scattering 

We now proceed to reconstruct the potential u, and ultimately the solution x of 
We rewrite equation p.7p as 

T (^7^-5 z) Z'^' — zz\^ 

7-\ X(jn; z) = R{z)x{m; z) —— , (7.1) 

a(zj \p + izj 

where R{z) = is the reflection coefficient. This equation defines a jump condi- 
tion along 3z = between ^, which is meromorphic in the open half-plane 3z > 0, 
and X, which is analytic in the open half plane 5z < 0. Coupled with the knowledge 
of the behaviour of these functions for large |z|, namely 

- (l + o(-\] as|z|^+oo,3z>0 (7.2) 



a(z) \ \z 

X(m; z) - ^1 + O ^ as |z| -> +00, 5z < 0, (7.3) 

this becomes an example of the classical Ricmann-Hilbcrt problem. To solve this 
problem we consider the Cauchy integral defined along the real z-axis: 

I -TT- ^TTTT^^C with 3z > 0. (7.4) 

We first evaluate this integral by considering a semi-circular contour in the upper- 
half plane SyC > 0. Let r+ denote the arc C = Me*" with < 6* < tt and M > 1. 
Along this contour we have 

1 f T(m; C) 1 r iMe'" / ^ / 1 



2TTi Jr^ a{C)iC + z) 2T:i Jq MeiO + z \ \M 

1 
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Since a has a finite number of simple zeroes at z = Zk in 9z > and C ~^ ^ cannot 
vanish in this region, the residue theorem gives 

T{m;C){C^Zk) 1 



/ = hm 



C^.. a{C){C + z) 2 

efcX(m;zfc) / p - izk 
k=i ^'"•^ ^ ^^'^ 



fe=i 

N 



where 

Cfc = bk hm — — = , 7.6 

C-^zk a(C) a'[zk) 

and wc have used the fact that T(m; z^) ~ bkxijn\Zk) ( ^+|fr) ' ^'^^'^ by 
Theorem (|6.ip ek depends on n. Wc may also calculate the Cauchy integral (|7.4p 
by using the jump condition (|7.ip . 

1 r°°xKC),^^^ pi^(C)xKC) fp^r,c, (7.7) 



27riy_^ (C + z) 27riy_^ (C + z) Vp + < 

and considering a semi-circular contour in the lower half-plane 3C "Si 0. Let F 
denote the arc C = Me^^ with -tt < 6* < and A/ > 1. Then 



and since C + will vanish in < 0, the residue theorem gives 

^pxKCl,^_,.^ xKC)(C + .)^i 

27rii_^ (C + ^) C^-^ (C + ^) 2 

= -X{m; -z) + - 

= -X{m;z) -t- i. 

Thus (fTT)) gives 

2 27riJ_^ (C + z) Vp + </ 
and by combining this with (|7.5p to eliminate / we obtain 

^ (^ + ^aO + 27rzi_^ (C + z) + 

(7.8) 

For the following results we introduce the convenient notation 



A := 



p — iz 



^p + iz ^ 

Lemma 7.1. The function A™ obeys the orthonormality condition 

A"-'"u;(z)dz = <5„„, (7.9) 
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where w{z) = ^^^^ 



p- 

Proof. This follows directly from using the substitution ~ tan^^ (p) ' 
Proposition 7.1. x ^■s expressible as the discrete integral transform of 
x(m;z) = l+T[A'-n, 

^-1 + ^ E ^KOA^-", (7-10) 



r— — oo 



where the kernel K{m,r) is independent of z. 

Proof. Substituting (j7.10p into (|3.5p and using p ~ iz = x+T' P + "i-z ^ y+\ gives 



m+l 



2p i^(™ + 2,r + l)A'-"-(2p + M™+i) ^ /^(m + l,r + l)A''-"' 

r— — OO r— — CX3 

m+l m 

-(2p + ii™+i) A'(m+l,r)A'-™ + 2p ^ A'KrjA'^-" 

r— — cxD r— — oo 

+ lim [(2p/^(m + 2,r)-(2p + ?/,„+i)i^(TO + l,r))A'--"-i] =Au„+i. 

For such a K{m,r) to exist for r <m we therefore require that it satisfy 

2p [K{m + 2, r + 1) + K{m, r)] - {2p + u,„+i) [K{m + 1, r + 1) + K{m + 1, r)] = 0, 

(7.11) 

subject to the boundary conditions 

2pK{m + 2, m + 2) - {2p + u,n+i)K{m + 1, m + 1) = (7.12a) 

lim K{m,r) = 0. (7.12b) 



r— ^ — oo 



Note that we assume that K{m,r) -> fast enough so that X^^l-oo '')! 
exists, since this is sufficient for (|7.10p to exist along 3z = 0. We now show that 
a solution to (|7.1ip with boundary conditions (|7.12p exists and is unique for all 
r < m. By defining FQ{m) := K{m,m) the boundary condition (|7.12ap gives an 
inhomogeneous first-order difference equation for Fq, 

2pFa{m + 2) - (2p + Um+i)Fa{m + 1) = Um+i, 

which can be solved to give the one-parameter family of solutions 

m— 1 / 

FQ{m) = -\ + A^ W (l 



2 — — oo 



2^ 



The convergence of this product follows from 

i— — oo ?'— — OO \ i— — oo / 

The second boundary condition ()7.12bp then imposes = 1 and so Fq is deter- 
mined uniquely. Now set r = m in (|7.1ip and define Fi{m) = K(jn,m — 1). The 
equation then reads 

F^im + 2) - (1 + ^)i^i(m + 1) = (1 + ^)i^o(m + 1) - i^o(m), 
2p Zp 
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which is an inhomogeneous first-order difference equation for Fi in terms of Fq. Tliis 
can be solved using a summing factor and then the boundary condition (j7.12bp 
will determine Fi uniquely. Now set r = m — 1 in (I7.11|) and define F2{m) :~ 
K{m, m — 2). This gives an inhomogeneous first-order equation for F2 in terms of 
Fi , whiclr can be solved uniquely in the above fashion. Thus by induction wc have 
a unique solution K{m, m ~ a) for all to, where a is any non-negative integer. □ 

Proposition 7.2. For L < m, the kernel K{m, L) defined by (|7.10p can be ex- 
pressed as 

1 /■+°° d7 

K{m,L) = — T[A''-'"]A"-^ . (7.13) 

271 p + tz 



Proof. 



— / T[A'-'"]A™-^— ^ 
2nJ^^ p + iz 

1 r+°° f 2n \ dz 



J2 K{m,r)- / Xr-^^—dz 

TT p^ + Z^ 



r— — oo 



= ^ K{m,r)5rL 

r— — Qo 

=K{m,L). 



□ 



By inserting (|7.10p into (|7.8p and simplifying one then obtains 

2; -t- Zfe \p-tlZk/ 



IZ 

r= — oo fc=l 



.i"oo t^^{z + Zk){p-izk)\p + izk) 2mJ_^ C + z\p + iCJ 

. 1 f; ^,„.,., r ^fSiL (£^)'.c^o. (7.14, 

27r«^^_ J_oo (C + z)(p-iC) \P + iCJ 



r— — 00 



We now perform the integral in (|7.13p in order to recover K. To do this we multiply 
the equation by 2^{p+iz) ' '^here L < m, and as we have assumed 32 > we integrate 
in z from — cx) to -|-oo along a path P-f just above the real axis. 

Lemma 7.2. For L < to, z^ purely imaginary with \zi;\ > and 3^ = we /laue 
the following results: 

>■-■■ ..^^(i^y- (7.15) 



{z + Zk){p + iz) p-izk\p + izk 



(z-kC)(p + «^) (p-iC) Vp + ^C 
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Proof. This follows by considering the poles of the integrand and using the residue 
theorem appropriately. □ 

By Lemma (|7.2p . taking the aforementioned integral of equation (|7.14p gives a 
discrete Gcl'fand-Lcvitan integral equation for K: 

m 

K{m,L) + B{L)+ ^ K{m,r){B{r-m + L)+B{r-m + L-l))^Q, (7.17) 

r— — oo 

where 

Proposition 7.3. Suppose that for fixed m, \B{L)\ < |i?(m)| for L < m, and 
X]-cx) P-xists. Then the discrete Gelfand-Levitan equation (|7.17p has a unique 

solution K{m, L) for L < m. 

Proof. We first prove existence. If we express K as 

i^(m,L) = ^i7,(m,L), (7.19) 

then this will solve (|7.17p if 

m 

Hj+i = - Hj{m,r){B{r-m + L)+B{r-m + L-l)) 

r— — oo 

Ho = -B{L). 

For L < m the recursion relation can be upper-bounded by 
\H,+,{m,L)\<2 \B{r)\\H,{m,r)l 

r— — oo 

and so emulating the inductive process in Lemma (j4.3p one obtains 

2^'(Er=-ooi^wi)' 



Hj{m,L) < 



Thus (|7.19p converges absolutely and uniformly for L < m. If Ki and K2 both 
solve (|7.17p then f] := Ki — K2 satisfies 

m 

|/3(m,i)|< 2\B{r)mm,r)\. 

r— — QO 

If |/3(m,L)| < TV for L < m then by the above inductive argument we obtain 



\f^{m,L)\<N- 



for any j > 0, thus (3 — 0. □ 
The function B is dependent on the quantities 

S {i?(n; () for 3C = ; ; ek{n)} (7.20) 
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which collectively comprise the scattering data. The entire inverse scattering pro- 
cess can therefore be summarised as follows: Given a potential Um defined on 71 = 
that satisfies 2p + w,„ > for all m and 

+00 

j = -CX} 

one may construct the the Jost solutions ip{m; z) and 'ip{m] z) using equation p.Sp . 
This leads to the knowledge of the functions a and b and thus to the n-dependent 
scattering data through equations (|6.6p and (|6.7p . The kernel K{m, L) can then be 
calculated via equation (|7.17p . and the n-dependent potential reconstructed from 
the boundary condition (|7.12p : 

'1 + Kim + l,m + 1) 



2p 



(7.21) 



1 + K{Tn, m) 

The solution of (jl.ip can then be found by first writing Xm.n = P^tl + + C + 
f{m,n), where C is arbitrary and / — > as |7ti| +00, and then integrating 
/(m + 2, 7^) - f{m, n) = Um+i,n- This gives 

Xm,n = pm + qn + C + ^ Um-r,n- (7.22) 
r>l,7' odd 

8. One- and Two-Soliton Examples 
The potential Wm+i along ji — for the onc-soliton solution p.ip is 

^p"(p2 - 1) 

Um + l = —2k-, rrr^ — r—, r- (8.1) 

(ylp™+2 + l)(AP™ + l) ^ 

Solving p.Sp for 5(771; z) gives the Jost solutions as 

ip{m;z) _ fAp'^'aiz) + 1 



(p - iz)™ V ^P™ + 1 

^(m;z) /^p£+a(z) 

(p -f iz)™ V Ap^ + 1 

where 



(8.2a) 
(8.2b) 



z + ifc 

6(z) = 0. (8.3b) 
At zi = 7fc we have ■(/'(tti; zi) = Aip{m; zi) so 61(77; z) = Ap^^ and 

ei = 2ifc^p^. 
Since b{n; z) = we have i?(z) = 0, and 

2kA 
p + k' 

The discrete Gel'fand-Levitan equation (|7.17p for K{in, L) then becomes 
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It is natural to assume the form K{m,L) = Q{m)pp p^, yielding 

A(p-i - 1) l + ^p™-V^ 



1 + Ap^p^^ 



1 + Ap^p"^ 



Thus by ([T:^ 



2p 



(1 + Ap^p^^Y - (1 + Ap^"- V^)(l + Ap;»+V^) 

(i + ^p™-V?)(i + M"+V?) 
M"-V^(pp-i)' 



2p 



(i+^pir-v?)(i+Aprv?) 

2fc 2fc 



where we have used the fact that 2n = 2k ( 



pm + + C + 



Finally by (|7.22p we have 
2fc 



1 + Ap^^p^ 

Two-soliton solutions to (jl.ip are given in [3] and [13] , and are again derived by 
Backliind transformation and Cauchy matrix approach respectively. The solution 
along n = takes the form 

2(fci + fcs) + 2fci^ip^i + 2k2A2P^2 



Xmfi ^ pm + C + 



1 + A,p;\ + A2p;^, + (l^) A^ap-p- 



where 



p + kj 



q + kj 



^p — kj J ' ^ 
From this one may calculate Um+i, and then solve (j3.5l) to determine the Jost 
solutions as 

ipim;z) 



1 + (frfr) ^^p7i + (f+if) -^^^"2 + a{z)A,A,p^,p^, 



(p — izY 



1 + A,p;\ + A2p;i + (f^) AiA2p™iP™2 



(8.4a) 



V^(m;.) _ + (f+il) ^iPpI + (frfr) ^^PP2 + A,A,p^,p;^, 



{p + izY 



i + Aip™+A2p;'^ + (|^) A,A,p^,p;i 



(8.4b) 



where 



a(z) 



(z — iki)(z — 1/02) 
(z + ifci)(z + zfe) 



6(z) =0. 

Thus again this corresponds to a reflectionless potential, but we now have two 
discrete eigenvalues zi ~ iki and Z2 ~ ik2. At each of these we have 

ei(n) 2ifciAip;\ e2(n) = 2ifc2A2Pg2, 
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and the solution of the Gel'fand levitan is 
l+K{m, m) — 

L + AiPp^ Pgl+A2Pp2 Pq2+(^fe7+fe7j ^1^2Ppi Pp2 PqlPq2 



i + AlPpiPgi + A2Pp2Pq2 + [kT+k; ) ^i^'2PplPp2PqlPq2 

After simplifying the expression (|7.2ip for Um,7i and using (|7.22p one then obtains 

2(fci + ^2) + 2/c2Aip™iP^i + 2/ciA2P™2P^2 

Xm,n = pm + qn + C + 



/ N 2 ' 

i + AxPp^Pq^ + A2Pp2Pq2 + (^fcT+i; J ^1^2PpiP5iPp2P92 



which is (after redefining constants as necessary) the two-soliton solution given in 
[3] and [13]. 

9. Aribitrary Reflectionless Potential 

We now consider an arbitrary reflectionless potential, that is one that satisfies 
h{z) = on S5z = 00 We write the N discrete eigenvalues as Zj = ikj, and make 
the further assumption that < kj < q for all j . The Gel'fand-Levitan equation 
reads 

N . N m 

*■(•». ^ + E ^< + E Y:j<^<-. op^--' - 0^ (is.i) 

From equation ()5.13p wc have 

(2 p + ttr)</7^(r;z fc) 
2(P^ + zlY 



y ^JPpzl^^] >o, (9.2) 



^r— — oo 



and so we set —icj =: Cjp°p^^, where Cj and p° are constants satisfying Cjp'j > 0. 
To solve (|9.ip we assume the form 

^ 1 

K{m,L) = -J2Q,im)—-p:plpl, (9.3) 
,/=i P + 

which gives a system of N equations of the form 

Q,{m) + y Q^M = c,- (9.4) 



By defining 



Q(m)^ := [Qi(m),...,QAr(m)] (9.5) 

:= [ci,C2, ...,CAr] (9.6) 
PZP^ P'n c, 

, '1 1 (9-7) 



'^Onc can show that if I I4.26D holds, this implies that 
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we may express (|9.4p in matrix form as 

Q{mf{I + M) = c^. 

Since M may be written as the Cauchy matrix A,yj ~ ^ multiphed on the 
left by the diagonal matrix of elements p°p^^Pq^, and on the right by the diagonal 
matrix of elements Cj , the determinant formula for the Cauchy matrix A^^j gives 

1 / l'<j ^ ■' ^ 

Thus every term in the expansion for dct(/ + M) will be positive, and so / + M is 
invertiblc. Thus by (^3)) 

K{m + 1, m + 1) = -c^(/ + M)-i (p/ - k)-^x, (9.8) 

where ~ indicates the shift m i-> m + 1, is the diagonal matrix with entries 
f^vj — ^vjkj and r is the vector with entries 

rj^PjPp^Pq- (9.9) 

Theorem 9.1. Any solution Xm,n of p.ip that gives rise to a reflectionless poten- 
tial Um+i satisfying (j3.4p can be expressed as 

x„,^ri =pm + qn + C - [I + A/)"^r, (9.10) 

where M is defined by (|9.7p with positive constants Cj, C is arbitrary and r is 
defined by ([gj]) . 



Proof. We show that ([97TOl) satisfies (fr2T|) with A'(m + 1, m + 1) defined by ([9 
Using the identities 

(/ + M){pl + k) - {pi + k){I + M) = rc^ 
{pi - k){I + M)-{I + M){pl - k) = rc^, 

we have 

K{m + 1,171+ 1)m™+i,„ = -c^(/ + M)-^TC^{I + My\pl - kY^y 

+ c^{I + M)-^vc^{I + M)-\pI - K)-h 
= - c^{I + My^ [{pi - k){I + M)-{I + M){pl - n)] {I + My\pl - K)~^r 

+ c^{I + M)-^ [{I + M){pl + k) - {pi + k){I + M)] {I + M)-\pI ~ K)-^r 
= -c^{I + Aiy^r + 2pc^{I + My^{pl - ny'^r - c^{I + il7)"^?, 
where we have used the fact that {pi + n){pl — ny^r = r. Thus 

(1 + K{m + 1, m + l))u™+i,„ = -2pc^{I + My^pl - k)-^? 

+ 2pc^{I + My\pI - ny^r 

= 2p{K{in + 2, m + 2) - A'(m + 1, m + 1)) 

as required. □ 
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Remark 9.1. The A^-soliton solution given in [131 is 

Xm,n = A {pm + qn + C - {I + Aiy^r) 

+ B(-l)"+" {pm + qn + D- [I + Aiy^r) , 

where — = 1 and C and D are arbitrary constants. Thus the solution (|9.10p 
is the particular case of this solutions with A = 1 and B = 0. This is due to the 
fact that wc have assumed x ^ pm + qn + C as m ^ ±oo. This is not restrictive 
however as the full general solution can be obtained by noticing that if Xm,n solves 
(|l.ip then so does 

Wm,n ■= Axm.n + B (- 1)"""*"" (a;,„,„ + const.), 

provided that - B^ = 1. Thus we find that our arbitrary reflectionless solution 
(|9.10p agrees exactly with that given in [TJ] , the only difference being that here we 
have restricted the parameters such that (|3.4p holds. 

10. Conclusion 

In this paper we have developed an inverse scattering transform for the LKdV 
(jl.ip . Using the 3D consistency of the equation we generated a Lax pair following 
[12] and used this to determine the discrete Schrodinger equation ([3.5p and its direct 
scattering problem. This was solved in section ([3]), in which results were proved 
rigorously and precise estimates were obtained. The discrete "time" evolution of 
the scattering data was calculated in Section ^ using the second Lax equation, 
and in Section ([7]) the inverse problem, which was posed along the real axis of the 
spectral parameter, was solved. Rather than giving the solution in terms of the 
Jost solutions we emulated the continuous Riemann-Hilbert approach and derived 
a discrete Gcl'fand-Levitan equation, whose solution is related to the potential by 
(|7.2ip . The solution of the lattice equation (|l.ip is then given by (|7.22p . In Section 
([8]) the one-soliton and two-soliton solutions were given as examples and these 
were found to correspond to reflectionless potentials with one and two discrete 
eigenvalues respectively, as expected. The arbitrary reflectionless potential case 
with N discrete eigenvalues was then shown in Section (|9]) to correspond exactly 
to the A^-soliton solution given in |13| . except with certain parameter restrictions. 

The inverse scattering transform works for all potentials Um.n satisfying the 
summability condition ([4.15[) . which is the discrete analogue of the integrability 
condition in the continuous case, and the positivity condition p.4p . This positivity 
condition is sufficient to prove that the poles of the tranmission coefficient are 
simple, however it is not necessary as was shown in a counterexample. A more 
precise restriction on the potential is therefore desirable, and will be addressed in 
future investigations. 
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